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Abstract 

String equations of the p-t\i generalized Kontsevich model and the 
compactified c = 1 string theory are re-examined in the language of 
the Toda lattice hierarchy. As opposed to a hypothesis postulated in 
the literature, the generalized Kontsevich model at p = —1 does not 
coincide with the c = 1 string theory at self-dual radius. A broader 
family of solutions of the Toda lattice hierarchy including these models 
are constructed, and shown to satisfy generalized string equations. 
The status of a variety of c < 1 string models is discussed in this new 
framework. 
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1 Introduction 



The so called "string equations" play a key role in various applications of 
integrable hierarchies to low dimensional string theories. The most funda- 
mental integrable hierarchy in this context is the KP hierarchy that 



provides a universal framework for dealing with many KdV-type hierarchies. 
String equations for "(p, q) models" of two-dimensional quantum gravity can 
be treated in a unified manner in this language. In contrast, the status of 
the Toda lattice hierarchy which is another universal integrable hier- 
archy, had remained relatively obscure until rather recent years. The Toda 
lattice hierarchy was pointed out to be an integrable structure of the one- 
and multi- matrix models [|13|, but these matrix models (matrix integrals) 
were only considered as an intermediate step towards the continuous (double 
scaling) limit to two-dimensional quantum gravity. 

In the last few years, the Toda lattice hierarchy has come to be studied 
from renewed points of view, such as c = 1 strings 0, |23|, ^ |2^, two- 



dimensional topological strings 0, |T^, |T0|, 0, ||, the topological CP^ sigma 



model and its variations related to affine Coxeter groups |]TT], ^. As 
opposed to the {p, q) models in the KP hierarchy, these are related to string 
theories with a true continuous target space. Our goal in this paper is to 
elucidate the structure of those string equations, in particular, those of c = 1 
strings in a more general framework. 

It will be instructive to recall the relationship between the (p, q) models 
and the KP hierarchy. String equations of these models were first discovered 
in the form of the Douglas equation ||^ 

[P,Q] = 1, (1.1) 

where P and Q are ordinary differential operators of the form 



P = dl + a2dl-' + --- + a. 
Q = d^, + hdr' + --- + K. (1.2) 

These two operators were later found to be related to an extended Lax for- 
malism of the KP hierarchy. The extended Lax formalism, developed by 
Orlov and his coworkers [^, is based on the ordinary Lax operator L and 
a secondary Lax operator (Orlov-Shulman operator) M. They are pseudo- 
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differential operators of the form 

oo 
n=l 

oo 

M = ^nt„L"-i+x+^t;„L-'^-\ (1.3) 

n=2 n=l 

and obey tlie Lax equations 

|i.[B„,4 '^^IB,.,M] (1.4) 

and the canonical commutation relation 

[L,M] = 1. (1.5) 

The operators are given by 

Bn = (L")>o, (1.6) 

where "( )>o" denotes the projection onto the space spanned by nonnegative 
powers of dx] similarly, we shall use "( )<-i" for the projection onto the 
space spanned by negative powers of dx- To reproduce the Douglas equation, 
we define P and Q as 

P = LP, 

Q = --ML^-P + ^^^^1-^ + 1", (1.7) 
p 2p 

require the constraints 

(P)<_i=0, (Q)<-i = 0, (1.8) 
and restrict the range of the time variables as 

tp+q = tp+q+l = ■ ■ ■ = 0. (1-9) 

The Douglas equation follows automatically from the construction of P and 
Q — (|1.8| ) and (|1.9| ) ensure that P and Q are differential operators of the 
required form. Note, in particular, that ( |1.8| ) give constraints on the Lax 
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and Orlov-Shulman operators, the first one being a reduction condition to 
the p-th generalized KdV hierarchy. Thus the (p, q) model can be understood 
as a "constrained KP hierarchy" . From the standpoint of the KP hierarchy, 
therefore, it is these constraints on L and M rather than the Douglas equation 
itself that plays a more essential role. In our terminology, "string equations" 
mean such constraints as ( p..8|) . 

This correspondence with the KP hierarchy allows one to use many pow- 
erful tools developed for the study of the KP hierarchy, such as the Sato 



Grassmannian, the Hirota equations, Wi+oo algebras, etc. [^, [l5|, ^ |l2 
Since this raises many interesting mathematical issues, generalizations to 
multi-component KP hierarchies have been also attempted |38 |. 



The case of (p, 1) (or (!,(?)) models has been studied with particular 
interest in recent studies. This is the case where the models becomes "topo- 
logical", i.e., describes topological strings of Ak {k = p — 2) type p. Unlike 
the other (p, q) models, string equation for these models can be solved more 
explicitly in terms of a "matrix Airy function" (or "Kontsevich integral") [0], 
and by virtue of this matrix integral construction, it is rigorously proven (at 
least for the case of p = 2) that the corresponding r function is a generating 
function of intersection numbers on a moduli space of Riemann surfaces [21 . 
The (generalized) Kontsevich integral is an integral over the space of N x N 
Hermitian matrices of the form 

Z{A) = C{A) J dMexpTT{APM - -^MP+^), (1.10) 

where A is another N x N Hermitian matrix, and C{A) is a normalization 
factor that also plays an important role. Kontsevich's observation is that 
this integral has two-fold interpretations. The first interpretation, revealed 
by "fat graph" expansion, is that this is a generating function of intersection 
numbers. The second is that this is a r function of the KP hierarchy, 

Z{A)=T{t), (1.11) 

where A and t are connected by the so called "Miwa transformation" : 

1 ^ 
t„ = -TrA-'^ = ^A-". (1.12) 

Here Ai, . . . , A^r are eigenvalues of A. 
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As for the Toda lattice hierarchy, our knowledge on string equations is far 
more fragmental, but simultaneously suggests richer possibilities. The most 
fundamental and well understood cases are the one- and two-matrix models 
[p!^ , H, String equations of these two models are quite different. Let us 
specify this in more detail. As we shall review in the next section, the Lax 
formalism of the Toda lattice hierarchy uses two Lax operators L, L and two 
Orlov-Shulman operators M,M ^5,^. These operators are "difference" op- 



erators obeying a set of Lax equations and a twisted canonical commutation 
relations of the form 

[L,M]=L, [L,M] = L. (1.13) 

(M should not be confused with the matrix variable M in the Kontsevich 
integral.) String equations are formulated as algebraic relations between the 
two pairs {L,M) and (L, M). Roughly speaking, string equations of the 
one-matrix model are written 

L = L-\ ML-^ = -ML, (1.14) 

and those of the two-matrix model are given by 

L = -ML, ML-^ = -L-\ (1.15) 

The point is that both {L^ML'^) and {L-^,ML) are a canonical conjugate 
pair, and that string equations are canonical transformations between them. 
The difference between the one- and two-matrix models is that whereas string 
equations (|1.14|) of the one-matrix model are coordinate-to-coordinate and 



momentum-to-momentum relations, string equations ( |1.15D of the two-matrix 



model mix coordinate and momentum variables. Note that the latter is also 
a characteristic of Fourier transformations. This is actually related to the 
fact that the Kontsevich integral is essentially a matrix version of Fourier 
transformations p5[ . 

In fact, most examples of string equations in the Toda lattice hierarchy are 
variants of the above two. For instance, string equations of c = 1 strings 0, 
|23| , 27 1 and their topological versions [H, 16, 10, 3|, Q are of the two-matrix 



model type. The generalized Kontsevich models, too, may be considered as 
a solution of the Toda lattice hierarchy obeying string equations of the two- 
matrix model type |T^. Meanwhile, string equations of the topological CP^ 
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model and its variants O, 113, p[ are of the one-matrix type. Furthermore, 



the deformed c = 1 theory in the presence of black hole backgrounds [28 



are known to obey more involved string equations, though this case, too, is 
essentially of the two-matrix model type. 

In this paper, we are mostly concerned with string equations of the two- 
matrix model type in the above sense. We re-examine the generalized Kont- 
sevich models and the c = 1 string theory in detail, and present a broader 
family of solutions that includes these two examples as special cases. This 
will clarify the status of a variety of c < 1 string models. For instance, we 
shall show that the generalized Kontsevich model at p = —1 does not re- 
produce the c = 1 string theory, but is rather related to the Penner model 

. This poses a question on recent attempts in the literature that treat 
two-dimensional topological strings as "A^, strings at k = —3" . 

This paper is organized as follows. Section 2 is a brief review on necessary 
tools and results from the theory of the Toda lattice hierarchy. Sections 3 and 
4 deal with the generalized Kontsevich models and the c = 1 string theory. 
In Section 5, our new family of solutions and string equations are presented. 
Section 6 is devoted to conclusion and discussion. 



2 Preliminaries on Toda lattice hierarchy 

We first present necessary tools and results in an /i-independent form 
(i.e., letting ^ = 1), and show an /i-dependent formulation in the end of 
this section. Throughout this section, s denotes a discrete variable ("lattice 
coordinate") with values in Z, and t = (ti,t2, • • •) and i = (ti,t2, • • •) two 
sets of continuous variables that play the role of "time variables" in the 
Toda lattice hierarchy. 

2.1 Difference operators 

The Lax and Orlov-Shulman operators of the Toda lattice hierarchy are dif- 
ference operators of the form 

oo 
n=0 
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oo oo 

M = 5^ nt^.L'^ + s + ^ 

n=l n=l 



T — n, p'^^ -U n, p "+2 

L — UqC + Un+lG , 

n=0 

CXD OO 

M = -Y^ntnL-'^ + s+Y.VnL^, (2.1) 

?i=l n=l 

where e*^^" are shift operators that act on a function of s as e^^" f{s) = f{s + 
77,). The coefficients £t„ and Vn are functions of (t, i, q), Un = Unit, t, s), 

etc. These operators obey the twisted canonical commutation relations 

[L, M] = L, [Z, M] = L (2.2) 

and the Lax equations 

^=[i?n,M], |^ = [i?„,M], (2.3) 

where the Zakharov-Shabat operators Bn and -B„ are given by 

Bn = (L")>o, B„ = (L-")<o, (2.4) 
and ( )>o,<o denotes the projection 

(^ a„e"^0>o = E «"e"^% (^ «ne"^0<o = E (2-5) 

n n>0 n n<0 

We call ( |2.2| ) "twisted" because it is rather the "untwisted" operators ML~^ 
and ML~^ that give canonical conjugate variable of L and L: 

[L,ML"^] = 1, [L,ML"^] = 1. (2.6) 
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Another important set of difference operators are the so called "dressing 
operators" of the form 

oo 

W^l+J2 W"^' , Wn^Wn{t, t, S) , 

n=l 
oo 

W^Wo+Y. ^ne"^^ , Wn = Wr,{t, t, s) (2.7) 

n=l 

that "undress" the Lax and Orlov-Shulman operators as 

oo 

n=l 

oo 

L = We'^'W-\ M = W{s - J2 nine-''^^)W~\ (2.8) 

n=l 



This does not determine W and W uniquely, and one can select a suitable 
pair of W and W such that the following equations are satisfied: 

dtn_ Otn 

= {We^^^W-')>oW, = -{We-^'''W-')>oW. (2.9) 

These equations of flows in the space of dressing operators can be "lin- 
earized" as follows. A clue is the "operator ratio" 

U{t, t) - W{t, l)-^W{t, t) (2.10) 

of the dressing operators W — W{t,t) and W — W{t,t). One can indeed 
easily see that U — U{t,t) satisfies the "linear equations" 

so that the flows in the space of U operators are given by simple exponential 
operators: 

oo oo 

C/(i,i) =exp(^i„e"^^)c/(0,0)exp(-^i„e-"^»). (2.12) 

n=l n=l 
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Furthermore, the passage from {W, W) to U is reversible. Namely, given such 
an operator U{t,t), one can solve the above "factorization problem" ( p.lO| ) 
to obtain two operators W = W{t,t) and W = W(t,t), which then automat- 
ically satisfy ( p.9|) . This is a Toda lattice version of Mulase's factorization 
problem for the KP hierarchy |26|. Actually, this factorization problem is 
solved by reformulating the problem in the language of infinite matrices. 



2.2 Infinite matrices 

Note the following one-to-one correspondence between the sets of difference 
operators and of infinite (Z x Z) matrices: 

J2 an(s)e"^^ « diag[a„(0] A" (2.13) 

n n 

where 

diag[a„(i)] = (a„(2)%). A" = (5ij_„), (2.14) 

and the indices i (row) and j (column) run over Z. (This A should not be 
confused with the finite matrix A in the partition function of generalized 
Kontsevich models!) The projectors ( )>o,<o are replaced by the projectors 
onto upper and (strictly) lower triangular matrices: 

{A)>o={e{j-t)a,,), {A)^o={{l-e{j-t))a,,). (2.15) 

All the Lax, Orlov-Shulman, Zakharov-Shabat and dressing operators have 
their counterparts in infinite matrices, 

L ^ L, M ^ M, L^L, M ^ M, W ^W, W ^W, (2.16) 

and they obey the same equations as their counterparts. For instance, dress- 
ing relation ( |2.8| ) turns into a matrix relation of the form 

oo 

L = WAW-\ M = W(A + ^ nt„A")W-\ 

n=l 

oo 

L = WAW \ M = W(A - ^ nt„A"")W"\ (2.17) 

n=l 

Here A is the infinite matrix 

A = (iSij) (2.18) 
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that represents the multiphcation operator s. Similarly, the difference opera- 
tor U(t,i) has an infinite matrix counterpart \J(t,t), and its time evolutions 
are generated by exponential matrices: 



U(t,t) =exp(5]t„A")U(0,0)exp(-5]t„A-"). (2.19) 

n=l n=l 

Factorization relation ( |2.10| ) is now converted into a factorization problem 
of infinite matrices of the form 

U(t, t) = W{t,t}^^W{t,i). (2.20) 

This is an infinite matrix version of the Gauss decomposition, because W = 
W(t, t) and W = W(t, t) are lower and upper triangular matrices. As in the 
ordinary finite dimensional cases, this Gauss decomposition can be solved 
explicitly by (an infinite matrix version of) the Cramer formula [Q. In 
particular, the matrix elements Wn(t,i,s) and Wn(t,i,s) can be written as 
a quotient of two semi-infinite determinants. The denominators eventually 
turn out to give r functions of the Toda lattice hierarchy: 

r(t,t, s) = det(uij{t,t) {-oo < ij < s)). (2.21) 

This gives a Toda lattice version of a similar formula in the KP hierarchy 
[ pl[l . The matrix elements Uij{t,t) of \J{t,i) are connected with their "initial 
values" Uij{0, 0) by 



(0,0)S'„(-t), (2.22) 

m,n=0 

where S'„ are the fundamental Schur functions 

oo oo 

^5„(t)A" = exp(5]t„A"). (2.23) 

n=0 n=l 

Eventually, the infinite matrix (a GL{oo) element) U(0, 0) persists as arbi- 
trary constants in a general solution of the Toda lattice hierarchy. It is thus 
the GL(oo) group that plays the role of the Sato Grassmannian in the Toda 
lattice hierarchy. 
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2.3 Constraints 

"String equations" in our concern are derived from linear relations of the 
form 

A{A, A)U(0, 0) = U(0, 0)A{A, A), (2.24) 

among the matrix elements of U(0, 0). Here A{A., A) and A{A., A) are linear 
combinations of monomials of A, A and A~^: 

A{A,A) = 5]a^„A™A", 

m.n 

A{A,A) = ^a^^A'^A", (2.25) 

m,n 

m runing over nonnegative integers and n over all integers. A fundamental 
result is the following |3^, ||, [28[] : 



Theorem 1 The above linear relations of U(0, 0) is equivalent to the con- 
straint 

A(M,L) = A(M,L) (2.26) 

of the Lax and Orlov-Shulman operators, where 

A(M, L) = ^ a„„M"^L", A(M, L) = a™„M"^L". (2.27) 

m,n m,n 

These Constraints may be interpreted as a fixed point condition under 
Vri+oo symmetries of the Toda lattice hierarchy ^, ^ . These matrices 
are in one-to-one correspondence with difference operators, 

A(A,A) ^ A(s,e^O, v4(A, A) ^ v4(s,e^0, (2-28) 
and give a closed Lie algebra with the fundamental commutation relation 

[A, A] = A ^ [e^%s\ = (2.29) 
This is essentially a (centerless) W^i+oo algebra with generators 

) = (AA-^)'=A"+'= ^ (se-^0''e^"+''^^^ (2.30) 
It is this W^i+oo algebra that underlies the above constraints. 
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These constraints can be further converted into hnear constraints for the 
r functions of the form 

X^r(t, t, s) = X^r(t, t, s) + const. r(t,t, s), (2.31) 

where Xa = XA{t,s,dt) and X^(t, s,9f) are hnear differential operators in 
t and t that represent VTi+oo symmetries acting on the r functions, and 
"const. " is a constant that also depends on A and A. This constant may 
emerge due to a nonvanishing central charge in the W^i+oo algebra of symme- 
tries acting on the r functions. Fixing this constant requires a subtle trick 
[0, We shall not deal with this issue in this paper. 



2.4 r functions in Miwa variables 

A clue connecting generalized Kontsevich models and the KP hierarchy ||21], |1| 
is the Miwa variable representation of the KP r function. If the matrix 
U(0, 0) is upper or lower triangular, a similar representation of the Toda 
lattice r functions |jl9[ is also available as we show below. (In fact, such a 
representation persists without this condition, though the relation to U then 
becomes more complicated.) 

Theorem 2 (i) Suppose that U(0, 0) is lower triangular. Then by the Miwa 
transformation 

1 ^ 



^n = -EV" (2.32) 

oft with arbitrary parameters Xi, the r functions T{t,i, s) can be written 



s~N~l 



_ det(u,_j{0,t,Xi) {l<i,j <N)) 
rit,t,s)= / ,^ ^ . n (2.33) 

det(^Aj ^ {l<i,3<N)j i=-oo 



where Uj{0,t,X) are the following generating functions of matrix elements 
u,j{0,i) o/U(0,t).- 



oo 



uj{0, t, X)=J2 A"*"'uii(0, i), (2.34) 



(a) Suppose that U(0, 0) is upper triangular. Then by the Miwa transforma- 
tion 

1 ^ 



E/^r (2.35) 



n ■ , 

1=1 
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oft with arbitrary parameters fi, the r function can be written 

det(^._.(t,0,/i,) (l<^,j<iV)) ^-^-1 
r{t,t,s) = \ /sir. . . ^ 11 2.36 



i=—oo 



where Ui{t, 0, /i) are the following generating functions of the matrix elements 
M,,(t,0) o/U(t,0).- 

oo 

Uiit,0,fi)=J2ui.j{t,0)fi^. (2.37) 

j=i 

The infinite products of Uu and Ua in tlie above formulae are interpreted 
as follows. Note that the matrices U(0, i) in (i) and U(t, 0) in (ii) are still 
triangular. The Laurent series Uj{0,t,\) and ^(t, 0, u) thereby take such a 
form as 

oo oo 

Uj{0,t,X) =Y.\-'-^u,,{0,i), Ui{t,0,fx) = J2mAt,0)fi^. (2.38) 

i=j j=i 



As opposed to an ordinary setting |]19l, we do not assume that the leading 



coefficients are normalized to be 1. The infinite product then arises. This 
factor is essentially the same as encountered in the treatment of semi-infinite 
determinant (|2.21| ) p3[| , and can be interpreted as: 

„ ( const. Uoo...Unn {n > 0) 

H = < const. {n = -1) (2.39) 

[ const. /(m„+i,„+i . . _i) {n < -2) 

with an overall renormalization constant "const. " 



2.5 ?i-dependent formulation 

An /i-dependent formulation of Toda lattice hierarchy can be achieved by 
inserting h in front of all derivatives in the previous equations as: 

^-4-' W-W' P.40) 

dtn Otn Otn Otn 
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The discrete variable s now takes values in hZ. Accordingly, Lax and Orlov- 
Shulman operators are difference operators of the form 



n=0 

OO CX) 

n=l n=l 

OO 

Li — uoe + Un+ie^ , 

n=0 

OO OO 

M = -^nt„L-'* + s+^^;„L", (2.41) 

n=l n=l 

where the coefficients are functions of {h, t, t, s), and obey the twisted canon- 
ical commutation relations 

[L, M] = ni, [L, M] = hi, (2.42) 

and Lax equations 

dM dM 



dtn ' dt 



n 



dM - dM 

= M], = [Bn, M]. (2.43) 

Bn and Bn are defined in the the same way as in the case of ^ = 1. 

Note that such an /?.-dependent Toda lattice hierarchy emerges if one 
starts from an /i-indcpcndcnt formulation and rcscalcs variables as 

tn h-Hn, tn ^ h'Hn, S ^ k'^S (2.44) 

and 

M h'^M, M h'^M. (2.45) 

This is indeed the case for solutions that we shall construct in subsequent 
sections. The ^-dependent formulation, however, also admits solutions that 
cannot be obtained by the above rescaling of {t, i, s). 
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3 Generalized Kontsevich models 



Our first example is the generalized Kontsevich models that have been stud- 
ied in the context of the KP hierarchy. We introduce "negative times" i and 
reconsider these models in the language of the Toda lattice hierarchy. The 
idea is more or less parallel to the ITEP-Lebedev group [|l^], who however 
treated these models as solutions of a "forced" hierarchy on a semi-infinite 
lattice Z>o. We rather attempt to interpret these models as solutions of the 
hierarchy on the bi-infinite lattice Z. 

3.1 Partition function as KP r function 

In an /i-dependent formulation, the partition function of the p-th generalized 
Kontsevich model is given by 

Z{A) = C(A) J dM exp Ti (^Af M - -^A/P+^), (3.1) 

where the normalization constant C(A) is given by 

\p _ \P 

CiA) = const. Ht r 

X (det A)(P-i)/2 exp ^"^ Trf -J^Ap+A (3.2) 

+ 1 ^ 

By the standard method using the Harish-Chandra-Itzykson-Zuber formula 
p5|l , the partion function Z{A) can be rewritten as a quotient of two deter- 



minants, 

^ det(^„,(A.) (l<.,j<iV)) 



where 



c(A) = const. fi'^P-^'^/^ exp h-U^—XP+^). (3.4) 

+ 1 ^ 
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The right hand side of the definition of Uj{X) is an integral over the whole 
real axis. By the standard saddle point method, one can show that Uj{\) has 
symptotic expansion of the form 



WiW~E^""Vi (3.5) 
i=j 

as A ^ +00. It is these asymptotic series rather than the functions Uj{X) 
themselves that are eventually relevant for the interpretation of Z{A) as a 
KP r function. Namely, we insert this asymptotic expansion into the right 
hand side of ( |3.3|) and consider it as a function (or, rather, formal power 
series) of t by the Miwa transformation 

n ^ n 

^n = -E\"" = -TrA-^ (3.6) 



3.2 Extension to Toda lattice r function 

We now extend the above construction to Wj's with negative index j. Ap- 
parently, this appears problematical, because the integral (over the real axis) 
becomes singular at yU = 0; actually, what we need is asymptotic series rather 
the functions Uj{X) themselves. As we shall discuss later in a more general 
setting, the asymptotic expansion as A ^ +00 is determined only by a small 
neighborhood of the saddle point at = A. The other part of the path 
of integration may be slightly deformed to avoid the singularity at = 0, 
or even cut off! With this deformation or cut-off, the integral itself will be 
varied, but the difference is subdominant and does not affect the asymptotic 
series. Thus upon suitably modifying the definition of Mj(A), we can show 
that Uj{\) for j < 0, too, has asymptotic expansion of the form 

CXD 

«.(A)~EA"^"V,- (3.7) 

We can thus define Uij for all integers, and hence a Z x Z matrix U = U(0, 0). 
According to the general constrution in the last section, this determines a 
Toda lattice r function r(t, t, s). 
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Since U(0, 0) is lower triangular, the r function should have a finite de- 
terminant representation by the Miwa transformation 



h ^ h 

^n = -EV" = -TrA-^ (3.8) 



n ^ n 

1=1 



One can easily see that the generating functions Mj(0,t, A) is given by 

/I oo 
rf/x/i-^"^ exp h-^ (A> - (3.9) 



n=l 



with the same normalization factor c(A) and the same path of integration 
(deformed or cut off in the aforementioned sense) as in the definition of 
Mj(A) = Mj(0, 0, A). The r(t, t, s) can eventually be written 

det(^,-.,,,(0,t-A.) (l<^,j<iV)) 

T{t,t,s) = 1 — — — ^. (3.10) 

det(Af ^-^-^(l<.,j<iV); 



(Note that, along with t and t, the lattice coordinate s, too, has to be rescaled 
as s — i> h~^s in an ^-dependent formulation of the Toda lattice hierarchy, s 
now runs over the ^-spaced lattice hZ.) 

Formally, this r function can be derived from the following matrix integral 
that extends Z{A): 

Z{A, i, s) = C(A) J dMexp h'^ Ti^A^M - j^M^^^ 

oo 

-Y^t^M-^ + shgM). (3.11) 

n=l 

3.3 String equations 

It is now rather straightforward to derive string equations — we first derive 
linear relations among Uij = Uij{0, 0), then convert them into relations among 
{L,M,L,M). 

Linear relations can be derived from linear functional or differential re- 



lations among Uj's. This is just to apply well known calculations [|1^, ^ to 
Mj's with negative indices as follows. 
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First, by integrating by part, 
A%(A) 

= c(A) J dnh-^(e-K^{h-^X^ij)) X /x-^-^exp^-^(--^/x^'+') 

= c(A) / dn(n{j + + x exp;i-^(A> — ^--iJ" 

= h{j + l)uj+ii\)+uj.pi\). (3.12) 



Here we have implicitly assumed that the path of integration connects two 
points at infinity, so that no boundary term emerges. If the path has a finite 
endpoint, such as the semi-infinite line [e, oo), e > 0, the boundary term is 
subdominant compared to the main asymptotic series arising from the saddle 
point at iJ, — X. Since we shall eventually derive linear relations among Ui/s, 
such subdominant terms are negligible. 

Similarly (but without integration by part), 

= c(A) f diihX'-P-^(n-'-'exph"\X^ii 

J p(yx\ p + i 

These relations imply the linear relations 

1 p — ]^ 

Uij-i = H-- + ^^)ui-p,j + Ui+ij (3.14) 
P ^P 

among the coefficients of asymptotic expansion. In terms of U = U(0,0), 
they can be rewritten 

UA = (--AA-^' + ;i^^A-^' + A)U. (3.15) 
P 2p 
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Finally, we replace 



A — >L,L, hA — >M,M (3.16) 

(taking into account the rescaling due to the presence of h) and obtain the 
following result on string equations: 

Theorem 3 The Lax and Orlov-Shulman operators of this solution obey the 
string equations 

LP = ML-^ + P, 

L = --ML'P + h^^^^L-P + L. (3.17) 
p 2p 

3.4 Generalized Kontsevich model at p = —1 

Let us consider the case of p = — 1. The string equations then become 

L'^ = ML-^ + L-\ 

L = ML + hL + L. (3.18) 

As we shall argue in the next section, these string equations do not agree 
with string equations of two-dimensional (or c = 1) strings (compactified at 
self-dual radius). 

This is a puzzling consequence of our analysis. Some of recent studies on 
the two-dimensional topological string theory |^ are based on the hypothesis 
that the topological model at k = —3 (which corresponds to p = — 1) 
gives the two-dimensional string theory. Our result implies that such an 
identification is somewhat problematical. 

The model at p = —1, as Dijkgraaf et al. 0] remarked, is rather related 
to the Penner model. The partition function for p = —1 is given by 

» oo 

Z(A,t,0) = C(A) J dMexph-^ Ti (^A-^ M -log M -Y.tnM-''). (3.19) 



n=l 



If t = 0, this function factorizes into a power of det A and Z{1, 0, 0), and that 
the second factor Z{1, 0, 0) coincides with the partition function of the Penner 
model Curiously, however, Dijkgraaf et al. identified this matrix integral 



with the c = 1 string partition function. Actually, as we have mentioned 
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above, this matrix integral does not correspond to the c = 1 string partition 
function. 

What physical meaning do the "negative times" i possess? It seems likely 
that t are nothing else but the coupling constants of "anti-states" that Mon- 
tano and Rivlis postulate in their topological interpretation of Ward iden- 
tities for (1, q) models. This is in fact a main idea that lies in the heart of the 
work of the aforementioned studies |22] on the two-dimensional topological 
string theory. Although we consider their interpretation of the = — 1) 
model rather problematical, their work contains many intriguing ideas. We 
shall return to this point in the final section. 



4 Compactified c = 1 string theory 

Our second example is the compactified c = 1 string theory formulated in 
the language of the Toda lattice hierarchy 0, ^ |2^. This string theory 
possesses a discrete parameter /3 = 1,2,..., and P = I corresponds to the 
case at self-dual radius. 



4.1 Partition function of c = 1 strings 

The partition function of compactified c = 1 strings becomes the Toda lattice 
T function with a diagonal U = U(0, 0) matrix of the form [0, ^ 

= i-n:)^-^y^ ^In'^ .r ^- (4-1) 



In fact, this is a somewhat simplified version of a true c = 1 string parti- 
tion function ||^. The above definition is substantially the same as used by 
Nakatsu ||2^ and leads to the same string equations (except that fi in his 



definition is replaced by —h here). String equations of the original c = 1 



string partition function are presented in Ref. [2S] in a more general setting 
including black hole backgrounds. 

Since U = U(0, 0) is a diagonal matrix, the r functions have two Miwa 
variable representations with respect to t and i. In both representations, 
the T functions are written in terms of a finite determinant including the 
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generating functions ^^(0, i, A) and ■Ui{t, 0, fi). By the familiar integral repre- 
sentation of the Gamma function, one can easily obtain the following integral 
representation of these generating functions: 

X / d/x/x(-^-i)/2-^-^exp^-^((A//x)^ + /31og//-X: 

n=l 



i_ 1 

I h 2 



POO 

X / dAA('^-^)/'+' exp /i-i Ux/^i}^ -piogX + J2 tnX'') ■ (4.2) 

-''^ n=l 

4.2 Integral representation of Kontsevich type 

If (3 = 1 (i.e., in the case at self-dual radius), the r function has a matrix 
integral representation of Kontsevich type. For simplicity, let us consider the 
case of s = h{N — 1) (though a similar matrix integral representation persists 
in a general case) . By the Miwa transformation 

^n = --E/^r = --^M-, (4.3) 

the r function T{t, t, h,{N — 1)) can be written 

dei(uN-i-i{t,^,Uj) (1 <i,j < N)) 
r{tX HN - 1)) = const. \ V ^- -1 <- /vT^ " ^'''^ 

As in the case of generalized Kontsevich models, this can be converted into 
a matrix integral of the form 

/<x> 
dAexp;i-^Tr(AM-^ -logA + ^i^A") (4.5) 



n=l 



with a suitable normahzation factor C{M). 
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Although very similar, this matrix integral is substantially different from 
the generalized Kontsevich model at p = — 1. The latter rather corresponds 
to (3 = —1. (Mathematically, the present construction is also valid for 
/3 = — 1, —2, . . ..) If /? = — 1, one can indeed derive a matrix integral repre- 
sentation of the form 

/oo 
c/Mexp;i"^Tr(A-^M-logM- ^t^M"") (4.6) 

in the Miwa variables 



n=l 



h ^ h 

^n = -E\"" = -TrA-^ (4.7) 



(We have put s = just for simplicity. A similar expression persists for a 
general value of s.) This is nothing else but the generalized Kontsevich model 
at p = — 1 given by ( ^.191 ). 

Thus, as opposed to the claim of Dijkgraaf et al. [0], we conclude that the 
generalized Kontsevich model at p = —1 does not correspond to the c = 1 
string theory of /5 = 1, but falls into the case of the strange value /5 = — 1. We 
shall reconfirm this puzzling conclusion in the language of string equations 
below. 



4.3 String equations 

By the recursion relation r(x + 1) = xr(x) of the Gamma function, one can 
easily derive the following algebraic relations for U = U(0, 0): 

UA-^ = UlA + h^ + l)A-PV. (4.8) 

Here A and A are the infinite matrices introduced in Section 2. To obtain 
the associated string equations, we have only to resort to the substitution 
rule 

A^L,L, hA^M,M. (4.9) 

We can thus reproduce the following result of Nakatsu (who derives this 
result by a somewhat different method). 
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Theorem 4 The Lax and Orlov-Shulman operators of this solution obey the 
string equations 



= (--M-h^ + lW, 
^ P 2/5 ^ ' 

= + + (4-10) 

As promised, this resuh clearly shows that string equations (p.l8|) of the 
generalized Kontsevich model at p = —1 agree with none of the c = 1 string 
theory with f3 > 1. The only possible value is (3 = —1. 

One can rewrite these string equations into linear constraints (W^i+oo con- 
straints) on the T functions. To do this, just take the n-th power {n = 
1,2,.. .) of both hand sides of the above string equations as 

L"^ =(■■■)", L"^ =(■■■)", (4.11) 



and apply the relation between nonlinear constraints (|2.26| ) on (L, M, L, M) 



and linear constraints ( p.31| ) on the r function. If /5 = 1, 2, . . ., one will thus 



obtain a set of Wi+oo constraints of the form 
dT{t, t, s) 



dtnf3 
dT{t, t, s) 

dtnB 



Xp^ri{t,S,dt)T{t,t,s), 

X^,„(t,s,9t)r(t,t,s), (4.12) 



where Xn,i3 and X^,n are linear differential operators as in ( [^.31| ) that give a 



W^i+oo symmetry of the Toda lattice hierarchy. If /3 = 1, this is exactly the 
W^i+oo constraints that Dijkgraaf et al. present in the former half of their 
paper. (Thus, as far as we understand, they deal with two distinct cases in 
the same paper — the /3 = 1 case in the former half and the (3 = —1 case in 
the latter half!) Meanwhile, ii (3 = —1, —2, . . ., the VTi+oo constraints are of 
the form 

n\p\tn\p\T{t,i,s) = Xi3.a{i,s,dt)T{t,i,s), 
n\(3\in\(i\T{t,i,s) = Xi3.a{t,s,dt)T{t,i,s). (4.13) 
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As a final remark, wc would like to note that if h is replaced by —h in 
the definition of Uij, the string equations become 




(4.14) 



In several papers, string equations of c = 1 or two-dimensional topological 
strings are given in this form. 



5 Synthesis — generalized string equations 

We now present our new family of special solutions and associated string 
equations. These solutions have two discrete parameters {p,p). The gener- 
alized Kontsevich models and the compactified c = 1 string theory can be 
reproduced by letting these parameters to special values. In fact, there are 
two apparently different constructions starting from the generating functions 
Uj{X) and Ui{fJ^), respectively, which eventually lead to the same string equa- 
tions. Unfortunately, we have been unable to see if these two constructions 
give the same solution. We mostly present the first construction based on 
Uj{X), and just briefly mention the second one. 



5.1 Generating functions and string equations 

The generating functions Uj{X) of the new special solution are given by 



u 



(A) = c(A) /dAi/^(P-^)/2-^-^exp;i-VAV —fJ^"^^ 

J \ p -\- p . 



c(A) = const. X^-^^/^exph-\^^XP+P). (5.1) 

The case of p -|- p = is also included here by replacing 

^ A^'+^"^plogA (5.2) 



p + p 



in the limit as p — >■ —p. In fact, Uj{X) then becomes monomials of A, and 
reduce to the generating functions in the c = 1 string theory. The case of 
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{p,p) = (p, 1), meanwhile, is nothing but the p-th generahzed Kontsevich 
model. The integrands of the above integrals may have singularities at the 
origin, but this apparent difficulty can be remedied by the same trick as 
discussed in the previous cases. Furthermore, in much the same way, one 
can derive the following relations for general {p, p) : 

The solution of the Toda lattice hierarchy in question is determined by 
asymptotic expansion of Uj{X). In fact, whereas the case of p + p > is more 
or less parallel to generalized Kontsevich models, the case of p + p < is 
quite distinct — we have to consider asymptotic expansion as A — ^> +0 rather 
than A +oo. To see this, let us change the variable of integration from n 
to a new variable z as: 

/i = \z. (5.4) 

Then 

Uj{X) = const. A(P+^^)/2-i-i exp ^-i(-^AP+P") 

X / e^p(^fi-^XP+P) ( t_zP+p) . (5.5) 

J ^ p + p ^ 

Thus the saddle point analysis has to be done in a region where h^^X^^P — >■ 
+00. The most dominant contribution comes from the saddle point z = 1, 
which corresponds to /i = A in the original integral. Contributions from other 
saddle points (and from the endpoint of the path of integration, if the path 
is cut off in a neighborhood of the origin) are subdominant and negligible 
in asymptotic expansion of Uj{X). The prefactor on the right hand side of 
( |5.5|) cancels the leading quasi-classical contribution from the integral, so 
that Uj{X) has asymptotic expansion in negative integral powers of A^"*"^. In 
summary, we obtain the following result: 

Theorem 5 (i) If p + p > 0, Uj{X) has asymptotic expansion of the form 

oo 

UjW ~ E (A ^ +oo). (5.6) 

i=j 



n 



j + i , p 



p 



2p 



-p 



-^"4 - » , 

p oX 2p 



- 1 



-x-p + xnuj{x) 



(5.3) 
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In particular, U is lower triangular. 

(a) If p + p < 0, Uj{X) has asymptotic expansion of the form 



j 



U,{X) ^ Y: >^''''^^J ^ +0)- 



(5.7) 



oo 



In particular, U is upper triangular. 

(Hi) In both cases, Uij — ifi—j^O mod p + p. 

We shall examine this asymptotic expansion in more detail in the next 
subsection. Anyway, the infinite matrix U = U(0, 0) determines a solution 
of the Toda lattice hierarchy. By the triangular form of U, the r function is 
ensured to have a Miwa variable representation. One should however note 
that one cannot use Uj{\) in the case of p + p < 0. If p + p < (so that 
U is upper triangular), it is Ui{^) that emerge in the Miwa variable repre- 
sentation. Unfortunately, we have been unable to give a closed expression 
to Mi(/i) in the present setting. [This should not be confused with the "dual 
construction" discussed later, which starts from an explicit construction of 
Ui{ii). In that case, the corresponding Uj{\) remains unknown in turn.] This 
is an incomplete aspect of the present construction. 

The aforementioned relations among the generating functions imply the 
following linear relations among the coefficients Uij of the asymptotic expan- 
sion: 




(5.8) 



In terms of the matrix U = (uij), 




1 



(5.9) 



String equations can be readily derived from these relations: 
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Theorem 6 The Lax and Orlov-Shulman operators of this solution obey the 
string equations 

LP = -ML-p -h^^L-'p + LP, 

p 2p 

LP = --ML-P + hJ^^L~P + LP. (5.10) 
p 2p ^ ' 

The string equations of generahzed Kontsevich models and c = 1 strings 
can be reproduced by letting {p,p) = (p, 1) and {p,p) = — The case 
of {p,p) = (1,-1) is also essentially the same as the string equations of 
the two-matrix model ; the extra term L and L on the right hand side can 
be absorbed into a shift of s. (More precisely, it is a "forced Toda lattice 
hierarchy" confined to a semi-infinite lattice thar arises here.) The string 
equations for general values of (p, p), too, relates "coordinates" L and L to 
"momenta" ML~p /p and ML~p jp. In this respect, these string equations 
may be referred to as "two-matrix model type" , as mentioned in Introduction. 



5.2 Systematic method of asymptotic expansion 

Let us examine the coefficients uij of the asymptotic expansion in more detail. 
We would eventually like to evaluate them in an explicit form. Such an 
explicit formula will be useful in considering the relation between the present 
construction and its "dual" form discussed later. 

Since this issue can be treated in a more general context, we now consider 
an integral of the general form 

l{k) = J dzg{z)e-^f^''\ (5.11) 

where /c is a large positive number, f{z) and g{z) are holomorphic functions, 
and f{z) is supposed to have no saddle point other than a simple nondegener- 
ate saddle point ai z = zq along the path of integration. By the assumption, 
f{z) can be written 

fiz)=f, + Mz-Zoy + 0{iz-z,f), /i^O, (5.12) 

in a neighborhood of zo- Then I{k) should have asymptotic expansion of the 
following form as k +oo: 

oo 

I{k) ~ e-'^^o «nfc"""^'/'^ (5.13) 

n=0 
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The problem is how to evaluate these coefficients. From the field the- 
oretical standpoint, the most familiar method would be the expansion into 
"Feynman graphs". Namely, one separates f{z) into the quadratic part and 
the rest, as shown above, treating the latter as a "perturbation", and evalu- 
ate the final series of integrals as Gaussian integrals. This is actually a very 
inefficient way, because a„ is then given by a sum of various Feynman graphs, 
and evaluating the sum is a hard task. A more efficient method is to resort 
to the Laplace method. A modernized version of this method is presented 
by Berry and Howls [H (whose main concern rather consists in global issues 
like Stokes phenomena). This method yields a closed formula for a„, as we 
now briefiy present below. 

A clue of the latter method is to rewrite the integral I{k) by the new 
variable 

y = {f{z) - = fl'\z - z,) + 0{{z - ^o)'). (5.14) 

Note that it is just a neighborhood of zq that actually contributes to the 
above asymptotic expansion; "modulo subdominant terms", one can replace 
the full path of integration by a small piece in a neighborhood oi zq. ?/ is a 
local coordinate in such a neighborhood of Zq. The integral I{k) can now be 
written into a Gaussian integral: 

Uk) = 6-''^° I dy^g(z)e-^y' + subdom. (5.15) 
J dy 

The subdominant terms "subdom." are of course invisible in asymptotic 
expansion, and negligible in the following calculations. The Jacobian dz/dy 
and the function g{z) are both holomorphic functions of ?/ in a neighborhood 
of y = 0. Let us write the Taylor expansion of their product as: 

dz °° 

g{z) = Y^h^y-. (5.16) 

The asymptotic expansion of the above integral can be obtained by inserting 
this Taylor expansion, formally interchanging the order of summation and 
integration, and finally extending the range of integration to the full real 
axis. Final integrals are to be evaluated by the familiar formula 

-OO ^ J —oo 
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Thus ttn are given by 

a„ = r(n+^)/i2n- (5.18) 
Furthermore, the Taylor coefficients /i2n are given by an integral of the form 

If , dz giz) . 
''-=2si„*d^siTI (8.19) 

along a small contour around y = 0. This contour is mapped onto a contour 
around z — Zq on the z plane, because y — > z is a coordinate transformation. 
Thus the last integral can be rewritten into a contour integral on the z plane: 

h2n = ^ / dz- (5.20) 

In summary, we obtain the following formula for the coefficients a„: 

+ r n(-) 
an = — 4> dz -^^'^ (5.21) 

In the case of our Uj{X), f{z) and g{z) are given by 

f{z) = -^^z'^' - ^\ = z^-'-^l'^-^-^ . (5.22) 

The relevant saddle point is at z = 1; the others are p-th roots of unity and 
the origin, and all subdominant if we consider the asymptotics as yc^yp^^ — > 
+00. The above contour integrals are still hard to evaluate explicitly, but 
anyway, we thus obtain a closed formula for the coefficients u^y 

5.3 Dual construction 

A parallel construction is possible by starting from the generating functions 
UiU) = c(/x) /dAA^f-'^/'+^exp^-^AV- 



c(ii) = const. //(^-')/'exp;i-V—^/iP+^^). (5.23) 

\p-\- p J 
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Previous calculations on Uj{X) can be repeated in a fully parallel manner. 
The generating functions Ui{n) turn out to have asymptotic expansion of the 
form 

i 

~ H (/^ +00) (5.24) 

j=-oo 

if p + p > 0, and 

00 

Uiifi) ^Y.Uijf^^ (a^^+O) (5.25) 
j=i 

if p + p < 0. (The case of p + p = again reduces to the c = 1 string 
theory.) The coefficients Uij satisfy exactly the same linear relations as the 
coefficients Uij in the previous construction. In other words, this "dual" 
construction leads to a solution of the same string equations. 

We have been, however, unable to prove (or disprove) that these two 
constructions in fact give the same solution. An obstruction is the fact that 
the linear relations among Uij appear to allow more ambiguities than the 
overall rescaling Uij — *• const. Uij. Nevertheless, it seems likely that the above 
solutions are rather special and coincide. The contour integral representation 
in the previous subsection will be useful in pursuing this issue. 



5.4 Concluding remark of this section 

The contents of this section is inspired by the work of Kharchev and Mar- 



shakov |20] on the (p, q) duality in c < 1 strings. Our U is a matrix represen- 
tation of their Fourier-Laplace integral operator that gives the (p, q) duality. 
Namely, if Vpg denotes the point representing the (p, q) model in the Sato 
Grassmannian, U acts on the Sato Grassmannian and interchanges the (p, q) 
and (g,p) models: 

VVgp = Vpg. (5.26) 

Note that everything is formulated in the language of the KP hierarchy. 

Although closely related with their work, our usage of U is substantially 
different. We rather interpret U as a GL(oo) element that determines a solu- 
tion of the Toda lattice hierarchy. Only in the case of g = 1, their work and 
ours are directly related via generalized Kontsevich models. Recall, however, 
that even in that case, our interpretation of generalized Kontsevich models 
is slightly different from the ITEP-Lebedev group |TP[ . We have pointed out 
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that the generahzed Kontsevich models with "negative times" give a solutions 
on the full (bi-infinite) lattice rather than a half (semi-infinite) lattice. 

6 Conclusion and discussion 

We have constructed a new family of special solutions of the Toda lattice 
hierarchy, and derived string equations of these solutions. These solutions 
have two discrete parameters {p,p), and include already known solutions as 
follows: 

• = (p? 1)) P > — the p-th generaized model with "negative 
times" i and a "discrete time" s. 

• (PjP) = iP: > — the compactified c = 1 string theory. [3 = 1 
corresponds to the self-dual radius case. 

• {PiP) = 1) — related to the Penner model. This case may be 
interpreted as the generalized Kontsevich model at p = — 1, and also 
as the c = 1 string theory "at (3 = —1". 

In the course of reviewing known examples, we have also posed several 
questions on the relation ^ between the generalized Kontsevich models 
and the c = 1 string theory. In particular, as opposed to the hypothesis 
postulated in the literature ||2^, the naive extrapolation of the p-th gener- 
alized Kontsevich model to p = —1 does not give the c = 1 string theory 
at self-dual radius. The p = —1 limit is rather related to the Penner model 
[0. Actually, these three models are mutually distinct and obey different 
string equations. Thus, although there is a lot of evidence that the Ak+i 
model at A; = —3 is anyhow related to c = 1 (or two-dimensional) strings, 
the approach to c = 1 by the generalized Kontsevich models seems to be 
problematical. 

A natural extrapolation of the c = 1 string theory at self-dual radius will 
be given by the case of p = 1 and p < —1. This model is ensured to reduced 
to the c = 1 model by letting p — >■ —1. Furthermore, remarkably, this model 
has a matrix integral representation of Kontsevich type, as follows. The 
matrix U = U(0, 0) of this model is upper triangular; the r function has a 
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Miwa variable representation by the Miwa transformation 



h ^ h 

tn = — E/^r" = — TrM-" (6.1) 



n—[ n 



This Miwa variable representation can be rewritten into a matrix integral of 
the form 



T 



{t,t,HN-l)) 



C(M) 



/■ 1 oo 

J dA exp h-^ Tr (AAF - ^^^^^^ + E ^^A") , (6.2) 



where we have considered the case of s = ^(A^ — 1) for simplicity. (This 
expression can be readily extended to general s.) 

Unlike the ordinary generalized Kontsevich models, the Kontsevich po- 
tential A^"*"^/ (1+p) in this matrix integral is a negative power of A. Therefore 
the integral should be suitably interpreted (e.g., by suitable regularization 
to avoid essential singularities as detA —* 0). The r function itself, as a 
function (or formal power series) of {t,i), is however independent of such a 
detail of regularization, and only determined by the contribution from the 
saddle point A = M. The above matrix integral representation might lead to 
a topological interpretation of this model by means of "fat graph" expansion. 

Another important difference from the ordinary generalized Kontsevich 
models is that the string equations of this model never reduces to string equa- 
tions in the KP hierarchy. In the case of generalized Kontsevich models, the 
"negative times" are coupling constants of unphysical states ("anti-states" 
in the language of Montano and Rivlis |^), which decouple from the theory 



as t — 0, and the string equations turn into string equations in the KP hi- 
erarchy in this limit. This does not occur in the case of p = 1 and p < —1. 
In other words, the model persists to be "c = 1" in the whole space of the 
coupling constants. 

We, however, do not know what physical interpretation the string equa- 
tions actually have for general values of {p,p) including the above case of 
p = 1 and p < — 1. Do they all correspond to a physical model (hopefully, of 
strings) ? A first step towards such deeper understanding will be to show a 
connection with moduli spaces of Riemann surfaces. The works of Montano 



and Rivlis [24] and Lavi et al. |22| will provide useful ideas in that direction. 
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Note added 

After completed this paper, the author learned that the same issue as we 
have considered on c = 1 strings was studied by Imbimbo and Mukhi in a 
recent preprint fl^. They noted the inconsistency between the Kontsevich- 
type representation of Dijkgraaf et al. ^ and the VTi+oo constraints, and 
proposed an alternative matrix model. Our matrix integral representation 
of the c = 1 partition function in Section 4 seems to be essentially the same 
as this model. The author would like to thank Sunil Mukhi for pointing out 
the presence of their paper. 
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